HEAVY QUARKONIUM PRODUCTION WITH POLARIZED 
HADRONS AND PHOTONS* 



MICHAEL KLASEN and BERND A. KNIEHL 
II. Institut fur Theoretische Physik, Universitat Hamburg 
Lumper Chaussee 149, 22761 Hamburg, Germany 
E-mail: klasen@mail.desy.de 

We investigate the inclusive production of prompt J/ip mesons in po- 
larized hadron-hadron and photon-hadron collisions in the factorization 
formalism of NRQCD. Numerical results are presented for BNL RHIC- 
Spin and the approved SLAC fixed-target experiment E161 to assess the 
feasibility to access the spin-dependent parton distributions in the po- 
larized proton and photon. We point out that data on J/tp production 
taken by the PHENIX Collaboration in unpolarized proton-proton colli- 
sions at RHIC tend to favor the NRQCD factorization hypothesis, while 
they significantly overshoot the theoretical prediction of the CSM. 

1 Introduction 

Before the advent of the non-relativistic QCD (NRQCD) factorization formal- 
ism, heavy quarkonium production with polarized proton or photon beams was 
believed to provide reliable information on the spin-dependent gluon distribu- 
tions of the polarized proton or photon. At present, however, we are faced with 
the potential problem that NRQCD predictions at lowest order have a consid- 
erable normalization uncertainty due to the introduction of non-perturbative 
color-octet matrix elements, which were not present in the color-singlet model 
(CSM) and which also have to be extracted from experiment [1]. 

In order to clarify the question if heavy quarkonium production with po- 
larized proton or photon beams remains to be a useful probe of the polarized 
gluon densities, we investigate inclusive J/ip production in polarized pp and 
72? collisions at RHIC-Spin and E161. Other processes with the potential to 
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constrain the polarized gluon distributions include the production of low-mass 
lepton pairs [2] . Numerical results relevant for J/ip production in polarized 77 
collisions at DESY TESLA as well as a complete list of analytical results for 
all polarized partonic cross sections can be found in Ref. [3] . 

2 J/i/> production at BNL RHIC-Spin and SLAC E161 

At RHIC, proton beams with longitudinal polarization of approximately 70% 
collide with center-of-mass energy up to y/~S = 500 GeV and luminosity C — 
2 x 10 32 cm _2 s _1 . In E161, circularly polarized photons with energies between 
35 and 48 GeV will collide on a fixed target made of longitudinally polarized 
deuterium [4]. In Figs. 1 and 2, NRQCD and CSM predictions for the unpolar- 
ized cross sections da/dpx and da/dy or da/dz are displayed in the first panel, 
while those for the double longitudinal-spin asymmetry All are shown in the 
second and third panels. The px distributions are integrated over the intervals 
\y\ < 2.4 and 0.3 < z < 0.8, respectively, while the y or z distributions are 
integrated over all kinematically allowed values of px in excess of 15 and 1.5 
GeV. The shaded bands indicate the theoretical uncertainties in the NRQCD 
and CSM default predictions. As our default polarized parton densities, we 
employ the GRStV-STD set in the proton and the GRSi-MAX set in the pho- 
ton [5] (for a detailed discussion of our input parameters see [3]). We assume 
the ideal case of 100% beam polarization. Realistic polarization is accounted 
for by scaling All with [-P(p)] 2 and P(j)P(D), respectively. 

At RHIC-Spin, the differences in All for various parton densities are 
large against the combined theoretical uncertainties from other sources, so 
that sufficiently precise measurements will increase our knowledge on the spin- 
dependent parton structure of the polarized proton. The NRQCD and CSM 
Pr-dependences incidentally almost coincide, so that the polarized proton den- 
sities can be explored in a model-independent fashion. On the other hand, the 
NRQCD and CSM predictions for the y distribution exhibit strikingly different 
shapes in the forward and backward directions. 

At E161, the situation is less favorable since the theoretical uncertainties are 
larger, due to the low photon-nucleon center-of-mass energy >/S = m^^E^ + 
mAr) rj 9.2 GeV. This is especially the case for the NRQCD prediction shown 
in the second panel of Fig. 2. Nevertheless, it should be possible to discriminate 
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Figure 1: Cross sections da/dpr and da/dy (first panel) and asymmetries All 
(second and third panels) of pp — > J/tfj + X at RHIC-Spin with y/S = 200 GeV 
as functions of px (left) and y (right) in NRQCD and the CSM. 
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Figure 2: Same as in Fig. 1, but for ■yD -> J/iji + X at E161 with £ 7 = 45 GeV. 
Inelasticity z is used instead of rapidity y. In the first panel, the resolved- 
photon contributions are also shown separately. 
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p p -> JA|/ X at Vs~= 200 GeV 




Figure 3: PHENIX data for B(J/t/> -> ^ + p~)d 3 cj/dy d 2 p T at RHIC with = 
200 GeV and in the interval 1.2 < y < 2.2 as a function of pr are compared 
with NRQCD and CSM predictions. 

between the GRStV and GS sets. In the CSM, the resolving power of All used 
to be much better, as is evident from the third panel of Fig. 2. In this sense, the 
introduction of NRQCD, which was necessary to overcome phenomenological 
and conceptual problems of the CSM, led to some aggravation. 

As can be seen from the first panels in Figs. 1 and 2, the normalization of 
the unpolarized cross section is a distinctive discriminator between NRQCD 
and the CSM. In fact, data from the PHENIX Collaboration [6] at RHIC, with 
y/S = 200 GeV, tend to favor the NRQCD prediction compared to the CSM 
one (see Fig. 3). Since the J/tp mesons are tagged through their decays to 
pairs, the factor B(J/ip — > = (5.88 ± 0.10)% is included in the 

theoretical predictions. We observe that, for pt > 2 GeV, the data are nicely 
described by the NRQCD prediction, while they significantly overshoot the 
CSM one. In the bin 1 GeV < pt < 2 GeV, the comparison has to be taken 
with a grain of salt since the NRQCD prediction and the P-wave contribution 
to the CSM one suffer from infrared and collinear singularities at px — 0, which 
still feed into that bin as an artificial enhancement. 
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3 Summary 

For inclusive production of prompt J/ip mesons in polarized hadron-hadron 
and photon-hadron collisions at RHIC and E161, we found that the spread in 
the asymmetries for different parton densities in general considerably exceeds 
the combined theoretical uncertainties from other sources, which we estimated 
rather conservatively. Therefore, even within the NRQCD formalism these 
experiments have discriminative power w.r.t. the spin structure of the proton 
and photon. As a by-product, we found that PHENIX data of unpolarized 
hadroproduction of J/ip mesons at RHIC favor NRQCD as compared to the 
CSM. This is in line with previous findings in pp, *yp, and 77 collisions [1, 7]. 
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